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r air  wanner  and  reduces  the  updraft  velocity  slightly  and  moist  downdraft  cools 
the  cloud  air  and  reduces  both  the  updraft  and  the  cloud  depth,  while  the 
influence  of  pressure  perturbation  is  also  to  reduce  the  updraft  but  it  is 
effective  only  for  clouds  of  radii  larger  than  1000  m.  The  influence  of  density 
stratification  is  to  reduce  wc  slightly  at  low  levels  and  augment  wc  at  high 
levels  and  to  slightly  increase  the  cloud  depth.  The  calculated  precipitation 
rates  appear  to  be  close  to  the  representative  values  obtained  from  observations- 
One  most  interesting  result  given  by  the  model  is  that,  under  the  less  unstable 
middle  latitude  environment,  the  precipitation  rate  jumps  from  lower  than 
.7  cm/hr  for  cloud  radii  smaller  than  1000  m  to  more  than  5.4  cm/hr  for  cloud 
radii  greater  than  2000  m,  reflecting  the  inhibiting  effect  of  the  mid- 
tropospheric  stable  layer  on  the  penetrative  power  of  the  smaller  clouds.  - 
Another  interesting  result  revealed  by  the  model  is  that,  with  the 
pressure  perturbation  included,  the  maximum  updraft  velocity,  wmax,  achieves 
highest  value  when  the  cloud  radius  is  about  8  km. 

The  influences  of  cumulus  convection  on  the  average  thermodynamic 
properties  of  the  atmosphere  have  also  been  calculated  and  it  was  found  that  the 
main  heating  is  from  latent  heat  but  compressional  heating  of  environmental  air 
is  also  of  importance,  especially  close  to  the  top  of  the  cloud. 

An  equation  for  the  influence  of  vertical  mixing  of  the  horizontal 
momentum  by  cumulus  convection  on  the  change  of  horizontal  momentum  has  also 
been  obtained. 
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1. 


Introduction 


It  is  well  known  that  the  total  life  span  of  a  vertically  oriented 
thunderstorm  or  a  single  precipitating  cumulus  cell  is  only  about  30 
minutes  and  it  can  be.  divided  into  three  distinct  stages,  namely,  an 
updraft  dominated  developing  stage,  a  mature  or  rain  stage  with  both 
updraft  and  wet-downdraft,  and  a  downdraft  dominated  dissipating  stage. 

It  is  also  known  that  the  strong  downdraft  within  the  domain  of  the  cloud 
itself  in  these  systems  is  attributable  to  the  downward  drag  created  by 
the  falling  rain  drops.  On  the  other  hand,  the  well  organized,  squall 
line  type  severe  convective  storms  usually  last  for  many  hours  even 
though  the  life  span  of  a  single  cloud  cell  within  the  system  may  also 
be  relatively  short.  It  is  evident  that  a  very  complicated  model  is 
needed  if  one  aims  at  representing  the  various  processes  Involved  during 
the  development  of  these  systems  in  detail  through  numerical  modeling. 
However,  essential  understanding  of  the  mechanisms  of  these  convective 
systems  can  be  achieved  by  the  use  of  much  simplified  models,  especially 
when  our  purpose  is  limited  to  the  finding  of  the  average  behaviors  of 
these  systems  only. 

In  this  study  we  shall  at  first  determine  the  representative  flow 
fields  during  the  whole  life  spans  of  the  cumulus  clouds  from  a  simple 
cloud  model  and  then  calculate  the  average  heating  effects  produced  by 
them,  and  finally  parameterize  the  average  heating  effect  in  terms  of  the 
large-scale  flow  variables.  Our  approach  to  the  cumulus  parameterization 
problem  is  still  as  that  discussed  by  the  first  author  in  an  earlier 
paper  (Kuo,  1974),  namely,  to  use  a  representative  model  cloud  or  a 
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combination  of  clouds  rather  than  use  an  unspecified  cloud  ensemble.  This 
approach  has  also  been  followed  by  Anthes  (1977)  in  his  approach  to  the 
cumulus  parameterization  problem,  but  with  a  somewhat  different  cloud 
model. 

Specifically,  the  cloud  model  we  adopted  in  this  study  is  a  quasi¬ 
steady  and  quasi-one-d imens ional  model  which  is  obtained  by  integrating 
the  two-dimensional  equation  for  the  azimuthal  vorticity  over  the  radius 
of  the  cloud  updraft  and  hence  it  takes  into  consideration  the  influence 
of  the  pressure  perturbation  approximately.  In  addition,  we  also  take 
into  consideration  the  influence  of  a  nearly  saturated  strong  downdraft 
in  the  immediate  environment  on  the  updraft  in  one  type  of  cloud,  which 
is  known  to  occur  often  in  nature,  in  addition  to  that  with  a  quiescent 
environment  only. 

2.  The  cumulus  cloud  model  and  the  governing  equations 

As  lias  been  mentioned  in  the  introduction,  the  flow  field  under 
consideration  is  taken  as  axisymmetric  and  hen  e  the  radial  and  vertical 
equations  of  motion,  the  heat  energy  equation,  the  water  vapor  equation 
and  the  anelastic  mass  continuity  equation  can  be  written  in  the  following 
forms  in  cylindrical  coordinates  (r,z): 

o 

ut  +  uur  +  =  -(p'  /PQ)  r  +  VV1U  ’  (D 

0 ' 

t  v  2 

wt  +  »«r  +  wwz  =  -(p  /Po>z  +  -  qc  )  +  vV  V,  (2) 

O'  +  r-i (urO ' )  +  p-1(p  wO')  +  wO  =  W20'  +  l*C/c  ,(3) 

t  r  o  o  z  o  z  p 

q'  +  r  ^(urq')  +  p  ^(p  wq ' )  +  wq  =  vV'q'  -  C,  (41 

^  r  o  o  o  z 

r  1 (ur)  +  P_1 (p  w)  =  0 
r  o  o  ^ 


i 


(5) 


} 


where  the  subscripts  t,  r  and  z  denote  partial  differentiations.  Here 

u  and  w  are  the  radial  and  the  vertical  velocity  components,  p',  9' 

and  q'  are  the  departures  of  the  pressure  p,  the  potential  temperature 

0  and  the  water  vapor  mixing  ratio  q  from  their  values  p^,  0  ,  and  q^ 

in  the  undisturbed  environment,  n  is  the  undisturbed  density,  0  = 

o  v 

T(1  +  0.608q)n  is  the  potential  virtual  temperature,  C  Is  the  rate  of 

condensation,  L  is  the  latent  heat  of  condensation,  q  Is  the  liquid 

water  mixing  ratio,  v  is  the  eddy  viscosity  or  eddy  conduction  coeffi- 

cient,  v1,  is  the  one -dimensional  Laplacian  operator  in  cylindrical 

coordinates  and  n  is  the  ratio  of  the  potential  temperature  to  the 

temperature,  which  are  given  by 
a 


V2  -  -4  +  -  ■  V,  -  V '  -  ! 

„  2  r  3  r  1 

3r  r 


.1000. R/ cn 

"  =  <-— )  P  , 


p  in  mb . 


(6a.b) 


(6c) 


The  undisturbed  quantities  p  ,  0  ,  q  and  p  are  taken  as  functions  of  z 

^  o  o  o  o 

only,  and  v  is  taken  as  proportional  to  w  in  the  cloud.  The  definitions 
the  symbols  used  in  this  paper  are  also  given  in  the  appendix. 

The  condensation  rate  in  a  unit  mass  of  cloud  is  taken  as  given  by 
the  formula 


C  = 


(7a) 


of 


where  qg  is  the  saturation  mixing  ratio  and  is  equal  to  q  in  the  cloud. 

If  all  the  liquid  water  are  of  the  nature  of  small  cloud  drops  which  move 
with  the  air,  then,  the  total  water  nixing  ratio  (q?  +  q)  of  an  individual 
parcel  without  nixing  or  prec i pi  tat i on  will  remain  constant  and  theroloro 


(  is  also  given  by 
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For  convenience,  we  add  7tL/c  times  (4)  to  (3)  and  obtain  the  following 

P 

equation  for  the  perturbation  entropy  s'  of  the  moist  air: 

5  S 

s'  +  r  1(urs')  +  p  1(p  ws')  +  w  r-^- 

t  r  o  o  z  3  z 


where 


=  fV2s  '  +  (q  -q  )w  , 

c^T  6 


s'  =  i  (0 
0 


.  +  ffLq'  ■ 


s  =  —  (6  +  - 

°  e  °  c. 


3  S  3  s  gLirq 
_ O  _  _ o  o 

9  2  “  32  “  c2TO 

P 


(8a,b  ,c) 


and  0  is  the  average  value  of  0  and  hence  is  a  constant. 

o 

For  the  saturated  region  inside  the  cloud,  it  is  more  convenient  to 

measure  the  mixing  ratio  perturbation  from  the  saturation  mixing  ratio 
» 

q  =  q  (T  )  at  the  undisturbed  temperature  T  rather  than  from  the  undis- 
'os  ‘so  r  o 

turbed  value  q  in  the  environment.  Since  q  is  equal  to  the  saturated 
o 

mixing  ratio  q  (T)  at  the  temperature  T  (=  T  +  T')  in  the  cloud,  we  have 

q  =  qsO')  =  qos  +  q"  »  (9a) 

where 

3qs  ,  9q 

q"  =  T'  =  —  — -)  o'  .  (9b) 

q  3  T  > it  3T  ;T  v  ’ 

o  o 

1'hus,  on  setting  q'  =  q"  ■+  q  -  q  in  s  ’  and  writing 

os  o 

-  I  <«’  +  ^  9")  -  f  [I  -  -  tVwv* 

P  pop 

s  =  =  (0  +  — ^-q  )  =  =6  =  s  +  (q  -q  )  ,  (9d) 

os  o  o  c  os  &  oes  o  c  n  os  's 
”  P  P 

we  find  that  Kq .  (8)  becomes 
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where 


Se>t  +  r_1(ursPr  +  +  W 


V';2se  +  “I?  (qs"qos)w  ’ 
c  0  1 
P 


3S  ,  10  dq  ds  gLuq 

OS  1  r _ o  +  Ltt  SO  1  o  & 

3z  n  dz  c  dz  dz  2- 

0  p  c  TO 

P 


(10a) 


We  apply  3/3r  to  (2)  and  3/3z  to  (1)  and  then  take  their  difference 

to  obtain  the  following  equation  for  the  azimuthal  vorticity  n  =  w  -u  : 

0’ 

\  +  (un^r  +  (wn)z  "  vVln  =  8  !7<F~V'  ~  \)  .  Ol) 

The  influence  of  the  pressure  perturbation  on  the  flow  field  Is  included 
in  this  equation  through  the  -u ^  term  in  n. 

According  to  Eq.  (5),  the  velocities  u  and  w  and  the  vorticity  n  can 
be  expressed  in  terms  of  the  stream  function  \p ,  viz.. 


u  =  -iji  +0ijj,  w  =  i|»  +  r  i>  , 

r 


(I2a,b) 


q  =  w  -u  =  V  ,4’  -ail i  +  ip 

r  z  1  Yz  zz  ’ 


( I  2c ) 


-1  3Po 

where  o  =  -p  is  the  density  stratification  factor.  Thus  the  dynamics 

o  3  z 

of  the  axisvmmetric  cloud  can  be  taken  as  governed  by  the  two  equations 
(10)  and  (11)  in  terms  of  the  two  dependent  variables  ip  and  s [  since  q^ 
and  hence  q"  of  (9b)  is  determinable  from  0. 
a.  Equations  for  the  cloud  updraft  region. 

Instead  of  solving  the  equations  (10)  and  (11)  in  the  two-dimensional 
space  (r,z),  in  this  work  we  shall  focus  our  attention  on  the  cloud  updraft 
by  reducing  these  equations  to  one-dimensional  forms  through  integration 


from  the  center  of  the  cloud  r=0  to  the  edge  of  the  cloud  r=R  under  the 
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assumption  that  we  have 

O'(R)  =  a0\  q'(R)  =  be'(R)  =  ab0',  (13a,  b) 

u(R)  =  -  o^w  =  -  Rw  ,  (13c) 

rR  > 

|  vrV  s^  dr  =  -  a"ws 1  ,  (13d) 

Jo  L 

where  a,  b,  at  ^  and  n"  are  constants  and  O',  w  and  s^  now  represent  their 
average  values  inside  the  cloud.  Here  a>  0  implies  entrainment  of 
relatively  warm  air,  and  b  >  0  implies  entrainment  of  air  which  is  more 
moist  while  h  <  0  is  air  drier  than  the  undisturbed  environment. 

In  this  study  we  used  a  =  -0.1,  b  =  0.002  for  the  moister  downdraft  and 
a  =  0.4,  b  =  -c  /Ln  for  the  drier  downdraft.  From  the  relations  (8a), 

(9c)  and  (13a,b)  we  also  find 


s' (Hi  =  n's' 
e 


s'  (R)  =  «'s'  -  (q  -q  ) 

e  e  cn  ns  o 

PU 


(13e, f ) 


whore 


3  % 

a'  =  a(c  +  bLn)/(c  +  L  y<=--)T  ).  (13g) 

P  P  3T 

We  now  multiply  (10)  by  p^rdr  and  integrate  over  r  from  0  to  R. 
After  making  use  of  these  simplifying  relations  at  the  cloud  boundary  we 


then  find  that  the  net  upward  transport  of  entropy  H(  =  p^wR^s^)  inside 


the  cloud  is  given  hv  the  following  equation: 


(H/w)  +  +  {K  -  ^J)H  = 


a*. 


;n 


ilz 


dS 

=  -  M  (T—  +  —  [  U  (q  -q  )  -  \  (q  -q  )]}  (14) 

dz  c  0  os  'o  c  T  .s  4is  ' 

P  P 

where  M(=  p  wR*^)  is  the  upward  mass  transport  inside  the  cloud  and  a*  = 
o 

2u^u',  where  o'  is  given  by  (13g). 

In  order  to  be  able  to  reduce  the  vorticitv  equation  (11)  to  one¬ 
dimensional  form  in  a  simple  manner  through  integration  over  r,  we  make 


k 


***** 


; 


use  of  Che  simplifying  assumption  that  the  radial  variation  of  ij/  is  oscil¬ 
latory  and  is  represented  by  a  function  of  r  only.  For  simplicity  of 
representation,  we  shall  follow  Holton  (1971)  by  equating  this  function 
with  the  first  order  llessel  function  J^Car/R)  with  a  -  2. AOS  Instead  of 
the  more  realistic  radial  function  given  by  the  linearized  equation 


(see  Kuo,  1965b).  We  then  have 


3_J|>  ,  1 

2  r  3  r 


(  1  5a ) 


Since  =  ---,  this  relation  also  implies 

7 

3w 

=  T  3~r  ' 

a 

Thus,  on  making  use  of  the  relation  (15b)  in  the  -u  term  of  n  we  then 


(15b) 


f  ind 


.  Ii  - 


2  \  9 

ct*"  0 


3)i  3-w  . 


(15c) 


Oz  3  r 

Further,  under  the  assumption  that  the  radial  and  the  vertical  variations 
of  w  are  represented  by  two  separate  functions  we  also  find  that  w t)  is 
given  by 


wtl  =  '2  ^--(wlw  -  — =•  (w  -  ow  )]} 
3r  2  zz  z 

a 

In  addition,  we  assume  that  we  have 


(13d) 


w(R)  =  0,  n(R)  =  -Pjh, 


vV  ndr  = 


« 

r  w  >1 
R  o 


(  1  6a  ,  b  ,  c  ) 


where  w  is  the  value  of  w  at  the  axis  and  q  is  the  mean  value  of  n 
o 


inside  R  and  is  given  bv 

R 


n 


1 

R 


2 

2  3  w 


3  w 


1  1  r  R  ,  o  o.i 

ndr  =  -  R  [wo  "  ~ 2  (7~2“  3T")] 


Ubd) 


JO  a  3  Z 

Observe  that  (16c)  is  of  the  same  form  as  (13d) 


Further,  we  assume  that  the  average  values  of  w,  and  q„  are  equal  to 
half  of  their  values  at  the  center,  viz.. 


Wo  =  2wc’  0v  (p)  =  2K'  Vo)  =  2q< 


(  1  6e ) 
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where  is  the  mean  updraft  velocity.  Then,  on  integrating  Eq  (11)  from 
r=0  to  r-K  and  making  use  of  the  relations  given  above  and  setting 

2  dw 

y  =  wc[wc  -  \  (— /  -o~)]  ,  (I6f) 

u  dz 

we  then  find  the  following  equation  for  y: 

t 

|t(y/wc)  +  +  ^)y  =  g[  (1-0. 5a)  f-  -  (qx  0.5qa  (R)J  ,  (17) 

R  L 

where  uj  =  Ujf  .  Thus  Che  nonsteady  development  of  our  quasi-one-dimen- 
sional  cloud  is  governed  by  the  equations  (14)  and  (17). 

We  point  out  that  in  this  formulation  of  the  problem,  the  flow 
conditions  in  the  environment  are  not  calculated  from  the  dynamic 
equations  but  their  influences  on  the  flow  conditions  in  the  updraft 
region  are  included  through  the  heat,  moisture  and  momentum  fluxes 
across  the  cloud  boundary  represented  by  (13a-d)  and  I6a-c)  and  through 
the  second  term  of  n  in  (15c)  which  represents  the  pressure  effect, 
b.  The  quasi -steady  equations. 

In  this  study  we  shall  limit  ourselves  to  the  quasi-steady  develop¬ 
ment  of  the  cumulus  convection  with  the  removal  of  rain  water  by  precipi¬ 
tation  included  but  without  giving  rise  to  a  downward  drag  on  the  updraft. 

This  can  be  taken  as  possible  either  by  assuming  that  the  ascending 
current  is  along  a  slanted  axis  or  assuming  rain  water  is  being  removed 
from  the  updraft  by  some  other  process.  Then  the  development  of  the  cloud 

can  be  perceived  as  a  gradual  extension  of  the  cloud  top  thro  eh  the 

successive  ascension  of  similar  buoyant  elements  originating  from  the 
cloud  base  until  the  equilibrium  level  is  reached,  while  the  decay  process 
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can  be  taken  us  accomplished  statically.  This  kind  of  development  can 
be  represented  by  the  steady  versions  of  equations  (14)  and  (17), 
namely,  with  the  time  rate  of  change  terms  omitted.  For  this  case  it  is 
more  convenient  to  solve  for  instead  of  the  net  upward  entropv  transport 
H  =  Msg-  Thus,  on  expanding  dfl/dz  in  the  steady  version  of  luj  (14)  and 
dividing  by  M  we  then  find  that  1L  reduces  to 

df> 

;)«’  =  - 

ir 


ds' 
e 

dz 


e  .  /  «*  .  a",  ,  _  os  tar  .  .  ,  g  ,  .  , 

v  +  mi  ~  p  +  ?)s.,  -  -  j',;  -  --  in(q  -q  )  -  (q  -q  ). 

z  K  o/  e  c  U  os  o  <•  I  s  os 

p 


p  o 


.-1 


where  p  -  M  dM/dz  =  2a  ^/R  Is  tlie  mass  entrainment  coefficient  and 
dS^/dz  Is  given  by  (10a).  When  expressed  in  terms  of  the  total  equivalent 
potential  temperature  0^(=  0  +  Lnq^/c  )  of  the  saturated  cloud  air  and 
when  bo tir  «*  and  a"  are  zero,  i.o.,  when  botlr  O'  and  q'  vanish  at  the 
edge  of  the  cloud  and  when  viscous  diffusion  is  neglected,  this  equation 
reduces  to  the  following  simple  form  : 
do 


.  -  +>*(0, 

dz  e 


«  )  =  ^  q 

oe  2„,  s 

e.  1 


(I  Ha) 


It  can  readily  be  shown  that,  when  expressed  in  terms  of  the  temperature 
T,  this  equation  takes  the  following  form: 

( 1  v  +  Lr-S-)  A'  +  t'[T-T  +  L-  (q  -q  )] 

dT  =  _  *  R  CpT _ °  Jj,  s  ° 


dz 


a + "-•‘“i 

c  p  dl 
P 


(18b) 


where  R  is  the  gas  constant,  is  the  virtual  temperature,  and  e  is  the 
saturation  vapor  pressure  at  the  temperature  T  of  tlie  cloud.  This  expres¬ 
sion  of  the  lapse  rate  in  tlie  cloud  is  roughly  equivalent  to  the  expression 
given  by  Stommei.  (1947)  but  is  slightly  more  general. 


(18) 


coefficient  p  as  known,  we  shall  use  this  equation  instead  of  (18a)  or 
(18b).  According  to  the  relations  (9a),  (9c)  and  (9b),  we  have 


1  (  \  ^_qs,  ®' 

0  (qs  q0s  3T  T  g 
o  D 


s'  3q 

— - _  (_^) 

3qo  -|  V3T  't  * 


r  l  3qs  i 

i+“. — —ii 
c  3  T  ;  T 

L  p  oj 


Therefore  (18)  can  also  be  written  as 


ds'  *  it  dS 

e  ,  ,  a*  a  \  ,  os  Lirp  .  . 

j-  +  (M  -  5-  +  ~ )«I  =  -  aZ - 


R  ’  2  e  dz  c  §  Mos  Mo' 

K  p 


where 


3cI, 


39e 


,«  =  u  _  _EL.  — 8/d  +  L__l£ 

P  2_,  3 T  c  3T 

c  I  p 

P  O 


>  , 


(19) 


(19a) 


dS  .  d0  T  dq 

os  _  1  r _ o  L_  ^os, 

dz  0  *dz  c  dz  * 

P 


(19b) 


Here  3q  /3T  is  for  T=T  .  Since  all  the  terms  on  the  right  hand  side  of 
so 

(19)  are  functions  T  and  p  ,  all  of  them  are  known  functions  of  z; 

o  o 

therefore  Eq  (19)  can  be  integrated  directly  provided  p’  is  taken  as  known 

For  q  and  3q  /3T,  we  make  use  of  the  following  relations: 

‘os  s 

0.622e  « 


qs  ~  P~e 


e^(T)  =  6.11  x  exp  {25.22  (1 


273K))  (273K)5.31 


rims,  from  these  relations  we  find 


3q  n  he 

s  _  0.622p  s 


3  T 


(P-es) 


2  dT 


(19c) 

(19d) 


(19e) 


The  variation  of  the  liquid  water  mixing  ratio  with  height  due  to 


condensation  is  calculated  from  a  similar  equation  In  accordance  with  (7a), 


viz.  , 


dq 


Jlcond. 


QZ 


dz 


U  {qs  -  q*  +  (l+a2)qjlcond>} 


(20) 


where  ^con(j  refers  to  the  total  liquid  water  mixing  ratio  created  by 
condensation  and  q*  is  the  mixing  ratio  at  r=R*  just  outside  the  updraft 
and  a^  is  another  constant  factor  for  allowing  the  effect  of  evaporation 
at  the  cloud  boundary.  The  formula  for  calculating  the  rate  of  precipi¬ 
tation  from  the  various  layers  will  be  given  later  in  section  4.  Thus 
our  system  is  represented  by  (19)  and  the  steady  version  of  (17),  viz., 

S  +  (iT  +  ^I)y  =  8  {(1-°*5a)  5-  -  q£  +  0.5  q*(R)}  •  (2]) 


C.  Two  types  of  immediate  environment. 

In  this  paper  we  shall  investigate  the  influences  of  two  different 
types  of  immediate  environment  on  the  updraft,  namely,  1)  the  quiescent 
far  away  environment  unaffected  by  convection  commonly  assumed  in  most 
one-dimensional  cloud  models  ;  2)  the  environment  characterized  by  a 

relatively  strong  downdraft  of  saturated  or  nearly  saturated  and  cooler 
air.  The  existence  of  such  saturated  downdraft  has  been  known  Jong  ago, 
and  their  cellular  patterns  have  been  revealed  more  clearly  by  more 
recent  Doppler  radar  observations,  see,  e.g.,  Ziegler  (1978).  On  tlie 
other  hand,  the  downdraft  in  the  far  environment  is  characterized  by  a 
much  slower  sinking  motion  and  a  temperature  higher  than  that  of  the 
undisturbed  state,  a  consequence  of  the  compressional  heating  in  the 
stably  stratified  environment. 
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The  colder  and  moist  downdraft  under  consideration  here  is  taken 
as  occurring  side  by  side  with  the  updraft  and  hence  its  nature  is 
different  from  the  downdraft  discussed  by  Srivastava  (1967)  and  Ogura 
and  Takahashi  (1971),  where  the  downdraft  is  directly  created  by  the 
faLling  precipitation  in  the  lower  part  of  the  original  updraft  and 
lienee  it  represents  the  decaying  otage  of  the  development  of  the  cumulus 
cloud.  One  way  of  creating  a  downdraft  side  by  side  with  the  updraft  is 
to  have  the  updraft  take  place  along  a  sloped  cold  front,  such  as  in  the 
squall  line  type  severe  storms.  Then  the  falling  precipitation  from  the 
updraft  can  create  a  downdraft  without  impeding  the  updraft.  Another 
way  is  through  the  outflow  of  the  raindrop  leaden  air  in  the  upper  layer. 
Since  our  simple  plume  model  is  not  able  to  describe  the  details  of  the 
development  of  such  downdrafts,  here  we  shall  simply  assume  that  such  a 
do'-ndraft  is  present  and  try  to  calculate  its  Influence  on  the  various 
properties  of  the  system. 

1,  Numerical  procedure 

As  has  been  pointed  out  before,  Eq,  (18)  can  t>e  solved  directly 
when  the  entrainment  coefficient  p  is  taken  as  known.  Therefore  the 
remaining  problem  is  to  solve  Eq .  (21)  by  making  use  of  the  values  of 
'V  obtained  from  (18).  Observe  that,  if  q^  is  also  a  known  function  of 
z ,  then  (21)  can  also  be  integrated  directly  as  a  first  order  equation 
for  y.  However,  as  described  in  section  4,  q^  is  a  function  of  wc . 
Therefore  (21)  must  be  solved  together  with  the  liquid  water  equation. 

The  procedure  we  adopted  in  solving  these  two  coupled  equations  is  an 
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Iterative  procedure,  which  is  first  to  solve  (21)  ievel  by  level  accord i nr 

to  the  traditional  plume  model,  that  is,  by  neglecting  the  pressure 

2 

perturbation  in  (21)  so  that  we  have  y  =  w  ,  and  using  w  (z  ,)  in  the 

c  c  n- 1 

liquid  water  equation  to  obtain  a  first  approximation  of  ,  and  use 

2 

it  in  (21)  to  obtain  y  =  w  for  the  next  level.  After  the  first 

c 

approximation  of  w  is  obtained  in  this  way  for  all  levels,  we  then 

recalculate  w  from  the  complete  relation  (16f)  which  includes  the 
c 


pressure  perturbation,  viz., 

.2 


n  ' 


d^w  dw 

c  c 

- - - ow  -  — 

c  ,  2  c  dz 

dz 


)  -  WC  =  -> 


(21a) 


This  equation  is  solved  by  Newton's  method  as  a  boundary  value  problem 

using  the  known  distribution  of  y,  including  possible  negative  values  of 

y  above  the  y=0  level.  Here  is  required  to  be  positive  everywhere 

except  at  the  top,  where  it  is  required  to  vanish,  while  at  the  cloud 

base  w  is  taken  as  known.  Eq .  (21)  is  then  solved  again  for  y,  using 
c 

the  w^  newly  obtained  from  (21a)  to  determine  a  new  set  of  q? .  The 
procedure  is  repeated  until  y  and  q^  are  no  longer  changing.  In  most 
of  the  calculations  presented  below  we  do  not  include  the  first  order 
differential  term  in  (21a)  which  contains  the  density  variation,  i .e. , 

,2 


.2 

R"  d  W 

"2  wc  T 
n  dz 


2  -  -  wc  =  -y 


(21b) 


However  as  revealed  in  a  later  section  it  is  of  minor  consequence.  In 
all  these  calculations  a  grid  step  equivalent  to  Ap  =  5  mb  is  utilized. 
The  location  of  the  upper  boundary  for  the  vertical  velocity  with 


pressure  perturbation  is  not  known  exactly  but  approximately  by  the 
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location  of  y=0,  hence  it  is  used  ns  a  first  guess  with  further  minor 
adjustments  made  during  the  calculations  to  give  a  best  fit,  i.e.  ,  to 
keep  positive  and  monotonically  decreasing  near  the  cloud  top.  As 
a  final  step  the  upper  15-30  per  cent  of  the  vertical  velocity  profile 
at  the  cloud  top  is  checked  by  recomputing  using  a  marching  procedure, 
making  use  of  existing  velocity  calculations  at  lower  levels.  Typically 
t  lie  only  consequence  Is  a  slight  adjustment  in  the  first  or  maybe  first 
two  noil  zero  values  of  next  to  the  cloud  top  while  at  all  other  points 
the  solutions  normally  agree  through  several  decimal  places.  This 

small  correction  is  made  necessary  since  the  location  where  w  =0  is 

c. 

usually  not  exactly  at  a  grid  point  but  it  can  be  found  from  the  marching 

solution  approximately  by  interpolating  linearly  between  the  last  positive 

and  first  negative  value  of  w  . 

c 

The  presence  of  a  pressure  perturbation  at  the  cloud  base  requires 

2  ? 

that  d  w  /dz*'  /  0  at  z=0.  An  estimate  of  its  value  can  be  obtained 
utilizing  a  backward  extrapolation  technique  built  around  Taylor  series 
expansion  and  a  one  sided  difference  approximation  of  (21),  i.e.. 


\y_  +  n v  ^  yB  / 


dz 

where  y  is  defined  by  (21b).  Assuming  that 


1  <JM  ,,  2 

M  dz  =  M  =  WcR  • 


the  approximation  can  be  shown  to  be 


(21*) 


2  2  2- 

d  w  (!+2ph)  w  .-(!-ph)w"  -hR  -  2w  ,w'  h+hR 

co  ci  co  o  cl  co 


dz 


2 


cl  __ 

2w  ,  h2  +  4  M  -M  ( 1  +ph) 
cl  2  1  0 

u 


(22) 


where  h  =  Az  and  the  subscripts  zero  and  one  refer  to  the  variables  being 

evaluated  at  z=0  and  z=h  respectively.  This  expression  quickly  converges 

2  2 

in  an  iterative  process  which  as  a  first  guess  assumes  d  w  /dz  =  0  at 

c 

z=0.  Thus  the  y  equation,  (21),  is  solved  with  an  initial  value  of  v 

'  n 

and  then  the  values  of  w^  are  found  from  Kq.  (21h).  Tills  process  is 

repeated  allowing  for  an  evaluation  of  F.q.  (22)  which  gives  a  correct  ion 

for  the  initial  value  of  y  until  the  change  in  y  at  z=0  is  less  than 

o 

2  -2 

.002  m  s  .  Then  the  procedure  is  as  described  above.  As  a  general  rule 

with  some  minor  exceptions,  it  is  found  that  clouds  of  small  radii,  say 

<  1000  m,  have  final  values  of  y  slightly  larger  than  w^  whereas 

o  ^  co 

2 

larger  radii  clouds  have  final  values  less  than  w  .  In  our  model  the 

co 

overall  consequence  of  these  changes  in  y^,  of  order  unity  or  less,  is 
very  small  hut  included  for  completion. 

Other  important  quantities  utilized  in  our  calculations  are  the 
fol lowing: 

entrainment  coefficient,  p  =  .183/Rj 

latent  heat  of  condensation,  L  =  597.0  -  ,555(T-273K)  CAL  gm  ; 

_2 

gravity,  g  =  9.814  ms  ; 
virtual  temperature,  T^  =  (1  +  .608q)T; 
and  finally  we  used  u"=0,  i.e.,  zero  viscous  diffusion. 
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4.  The  calculation  of  liquid  water. 

Liquid  water  exerts  a  negative  buoyancy  upon  the  updraft,  as  evident 

in  the  effective  buoyancy  term,  Eq .  (21),  so  its  concentration  must  be 

computed  simultaneously  with  the  vertical  velocity.  Liquid  water,  q^,  is 

really  the  sum  of  cloud  and  rain  water  where  cloud  water,  denoted  by  q^, 

is  composed  of  small  droplets  that  are  carried  by  the  updraft  and  rain 

water,  denoted  by  qf  ,  is  composed  of  the  larger  drops  that  fall  relative 

to  the  updraft.  Condensation  loads  the  cloud  via  cloud  water  whereas 

precipitation,  i.e.,  the  fallout  of  rainwater,  unloads  the  updraft. 

This  unloading  is  important  lest  the  :1'-ag  force  becomes  too  unrealistic. 

A  detailed  relationship  between  q  and  q  has  been  outlined  by  Kessler 

cw  rw 

(1969).  Here  a  simpler  format,  similar  to  that  used  by  Anthes  (1977) 
based  on  a  median  dropsize  is  followed. 

Cloud  water  is  known  to  increase  due  to  the  condensation  process 
and  to  decrease  because  of  the  autoconversion  process  which  creates  raindrops 
and  by  collection  of  cloud  droplets  by  raindrops.  Thus  in  symbolic  form. 


Aq  =  Aq  .  -  (Aq  )  „  -  (Aq  >  ... 

cw  c.ond  cw  auto  cw  coll 

The  autoconversion  process,  following  Kessler  (1965),  is  given  "by 


(23) 


(  q  ) 

Vw  auto 


KAt  (q  -a') 
Cw 


q  1-  a 
cw 


q  <  a 
cw 


(24) 


where  a'  =  0.5  gm  m  *  is  the  threshold  value  and  K  =  10  \s  ' .  The  collection 
term  is  expressed  as 


(Aq  )  -  5.26  x  10  *  q  q^'^’  At 

cw  co I  1  cw  rw 


(25) 


Rainwater  increases  at  the  expense  of  cloud  water  and  decreases  due 
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to  fallout,  thus  symbolically. 


Aq  =  (Aq  )  +  (Aq  )  .  .  -  (Aq  ),  (26) 

rw  cw  auto  cw  coll  ^rw  fallout 

It  will  be  assumed  that  the  fallout  rate  is  proportional  to  both 


q^  and  the  terminal  velocity  Vp,  so  that  we  have 


(Aqrw)  Fallout  XqrwVFAt_ 


(27) 


where  \  is  an  empirical  fallout  coefficient  whose  dimension  is  1/length 
and  hence  it  is  often  written  as  1/H  (Anthes  1977  ,  Krietzberg  and  Perkev , 
1976)  .  Here  both  and  q  are  expressed  in  gm  in  and  the  t  ime  step 
At  is  computed  from  Az/w^  for  w^  >  0.005  ms  \ 

Following  Simpson  and  Wiggert  (1969),  we  take  V{.  for  the  median 
size  rain  drops  as  given  by  the  following  formula  in  terms  of  q 

rw 


v  r  ,  0.L25  -1 

V  =  5 . 1  q  ms 

F  rw 


(28) 


-1 


The  rate  of  precipitation,  PR(gm  s  m  ")  per  square  meter  can  be 
computed  using 

N 

(29) 


PD  =  F.  At,  V„(Aq  )t.  / 1  , 

R  ,  i  F  ^rw  Fallout 
i  =  l 


where  At . 

l 

Ap  and 

t  = 


is  the  time  since  development  of  the  i 


th 


ln"er  of  thickness 


N 

>.  At  . 
i 


is  the  total  cloud  development  time.  It  is  eustomarv  with  plume  models 
to  equate  total  cloud  lifetime  and  development  time. 
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3.  Discussion  of  Results. 

In  this  study  the  characteristic  properties  of  cumulus  convection 
are  investigated  by  the  use  of  a  simple  quasi-Lagrangian  axi-symmetric 
steady  state  mode]  based  on  the  two-dimensional  vorticity  and  heat  equa¬ 
tions  which  include  the  effects  of  both  the  pressure  perturbation  and  the 
influence  of  downdrafts  in  the  environment.  These  equations  are  reduced 
to  one-dimensional  forms  by  integrating  them  from  the  center  of  the  plume 
outward  into  the  env i ronment .  The  entrainment  terms  of  Stommel  (1947) 
in  which  the  interaction  between  the  cloud  and  environment  is  represented 
by  an  expression  inversely  proportional  to  the  cloud  radius  are  recovered 
in  this  manner.  In  our  formulation  additional  terms  are  introduced  when 
nearby  environment;) 1  temperature  and  relative  humidity  profiles  differ 
from  the  undisturbed  state  far  removed  from  the  cloud.  These  relatively 
nearby  environmental  variations  owe  their  existence  to  the  updraft  and 
their  inclusion  into  the  model's  physics  is  distinctly  different  from 
tiie  lateral  mixing  utilized  by  Ogura  and  Takahashi  (1971).  The  influence  of 
the  compensating  nearby  downdrafts,  transport!  ,  s  of  lic  it  and  mass,  upon 
the  total  cumulus  convection  has  been  shown  to  be  of  major  importance  in 
several  observational  studies  (Johnson,  1976;  Nitta,  1978;  Zipser ,  1977, 
others).  According  to  Nitta  (1978),  during  the  active  stage  of  the  deep 
cloud  clusters  strong  downdrafts  exist  with  vertical  velocities  (amplitudes) 
about  40-50Z  of  the  adjacent  updrafts. 

Two  different  undisturbed  environmental  situations  are  investigated 
by  our  model,  namelv,  1)  a  potentially  more  unstable  environment  repre¬ 
sented  by  mean  Gulf  of  Mexico  hurricane  season  temperature  and  relative 


humidity  profiles  (Herbert  and  Jordan,  1959)  and  2),  for  later  considera¬ 
tion,  a  less  unstable  sounding  in  the  middle  latitude.  Pig.  ]  illustrates 
the  potential  (solid  line)  and  equivalent  potential  temperature  (dashed 

and  dotted  curves)  of  the  mean  tropical  sounding.  Two  0  ..  curves  are  given 

oh 

in  this  figure,  the  very  small  differences  between  them  are  due  to  the  fact 
that,  in  the  dashed  curve  the  latent  heat  of  vaporization  is  taken  as  con¬ 
stant,  namely,  L  =  bOO  cal  gm  \  whereas  in  the  dotted  curve  I,  -  597  -  .555l 
is  used,  where  t  is  in  °C,  and  the  latent  heat  of  freezing,  KO  cal  gm  is 
also  included  linearly  between  -10  and  -40°C.  Pig.  2  shows  the  vertical 
derivative  of  the  entropy  for  the  constant  L  case,  which  appears  as  a  lor 

function  in  the  perturbation  entropy  Eq.  (19).  A  positive  contribution  by 
this  term  is  indicated  below  595  mbs.  Unless  stated  otherwise  all  figures 
refer  to  calculations  using  this  mean  tropical  sounding. 

Here  we  shall  present  the  results  based  on  (21b)  first,  that  is, 
with  o=0.  The  cloud  temperature  excesses  obtained  from  our  plume 
model  for  a  cloud  with  R  =  1000  m  are  illustrated  in  Pig.  1.  Here  the 
dashed  curve  is  with  moist  downdraft  while  the  solid  curve  is  without 
moist  downdraft.  The  relative  humidity  of  the  downdraft  is  taken  as  95 
per  cent.  It  is  seen  that  the  downdraft  air  is  cooler  than  the  undis¬ 
turbed  air  represented  by  the  mean  sounding  and  as  a  consequence  the  over¬ 
all  depth  of  the  cloud  is  reduced  by  nearly  60  mb  in  pressure  measure. 

As  shown  in  Fig.  4,  the  correspond ing  reduction  of  the  vertical  velocity 
is  substantial.  Both  of  these  changes  are  the  results  of  the  reduction 
of  the  buoyancy,  which  itself  is  directly  attributable  to  the  fact  that, 
even  though  the  moist  downdraft  temp  •nature  T  is  lower  than  the  tem¬ 
perature  T  of  the  undisturbed  environment,  the  virtual  temperature 


perturh.it  ion  (T  -T  )  is  positive  on  account  of  the  higher  mixing  ratio 
Vdw  o 

of  the  moist  downdraft.  This  tendency  can  be  reversed  if  T.  becomes  much 

1  dw 

lower  than  T  .  Our  calculations  show  that,  for  the  same  temperature  and 
humidity  distributions,  the  influence  of  the  downdraft  on  the  buoyancy 
decreases  with  increasing  cloud  radius  and  vice-versa.  These  results 
ire  not  significantly  affected  by  the  presence  of  other  physical  pro¬ 
cesses  such  as  evaporation  at  the  cloud  boundary  or  inclusion  of  latent 
heat  of  fusion,  as  Illustrated  by  the  dotted  curve  in  Fig.  2  and  dashed 
curve  in  Fig.  1,  respectively.  The  latter  effect  is  introduced  by 
incorporating  the  latent  heat  of  freezing  linearly  between  0  and  -10°C. 

The  consequences  of  environmental  compressional  heating  upon  cloud 
P'lential  temperature  excesses  when  no  nearby  moist  downdraft  is  present 
are  displayed  in  Fig.  6  for  three  cloud  radii.  The  solid  lines  represent 
Llie  cloud  perturbation  in  the  absence  of  compressional  heating.  The 
overall  effect,  as  revealed  in  the  dashed  curves  (with  a=0.4),  is  to  warm 
the  updraft  especially  at  upper  levels.  The  vertical  velocity  however 
is  reduced  because  of  the  net  reduction  in  buoyancy  due  to  the  term  -a/2  in 
l.q.  (21)  that  more  than  compensates  for  the  enhanced  cloud  temperature 
excesses.  This  is  clearly  seen  in  Fig.  7  by  comparing  the  profiles  computed 
with  (dotted  curve)  and  without  (dashed  curve)  compressional  heating. 

For  a  small  radius  cloud,  R  =  500  m  (Fig.  7),  the  large  changes  in 
the  vertical  velocity  due  to  compressional  heating  dwarfs  the  influence 
due  to  pressure  perturbation.  The  latter  causes  a  very  minute  reduction 
as  indicated  in  Fig.  7  by  the  differences  between  the  solid  (no  pressure 
perturhat  ion)  and  dashed  curves.  However,  increasing,  the  cloud  radius  to 
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5000  m  substantially  enhances  tlu-  role  of  the  pressure  perturbation  so 

that  these  two  influences  become  nearly  equivalent  in  importance.  The 

effect  of  the  pressure  perturbation  as  revealed  in  this  figure  and  in 

others,  shown  and  not  shown,  is  to  reduce  the  magnitude  of  the  vertical 

velocity,  to  shift  the  maximum  vertical  velocity,  w  ,  downward  and  to 

max 

slightly  Increase  the  cloud's  depth,  as  can  be  seen  from  a  comparison 
between  the  solid  and  dashed  or  dashed-dot  curves  in  Pig.  H ,  while  the 
dotted  curve  illustrates  the  influences  of  compressionnl  heating  in  tin- 
environment.  The  distinction  between  the  dashed  and  dashed-dot  curves, 
both  of  which  contain  pressure  perturbation,  is  that  in  the  latter  the 
a  term  in  Eq .  (21a)  is  included.  The  influence  of  this  term  is  to 
further  reduce  the  velocity  by  a  small  amount  except  near  the  cloud  top 
where  w„  is  enhanced  by  it.  The  value  of  w  is  only  sl.ight.lv  decreased 
but  shifted  upward  in  this  deeper  cloud.  The  general  trend  with  or  without 
the  inclusion  of  this  first  order  term  is  much  the  same  hence  this  term 
is  ignored  in  most  of  oar  calculations.  These  conclusions  about  tin- 
pressure  perturbation  are  in  major  agreement  with  those  of  Ya"  (1179) 
obtained  using  analytical  models. 

Holton  (1974)  has  found  by  the  use  of  an  one-dimensional  Euler ian 
cloud  model  that  for  clouds  with  radii  greater  than  1  km,  the  pressure 
perturbation  substantially  reduces  the  growth  rate.  Within  the  frame¬ 
work  of  our  Lagrangian  model  it  is  not  possible  to  compute  the  growth 
rate.  On  the  other  hand,  another  very  Interesting  effect  of  the  pressure 
perturbation  is  revealed  by  our  solutions,  namely,  when  the  pressure 
perturbation  is  included,  achieves  its  highest  value  when  the  -loud  t.idius 


. •  t >> > ■  i L  K  km,  whereas  ill  tin-  absence  of  tin-  pressure  per t urb.'i t  i on  w 

max 

increases  monolon icallv  with  increasing  radius,  as  can  he  seen  from  the 
sol  id  and  the  dashed  curves  in  Fig.  9.  Our  results  also  show  that  small 
clouds  are  only  slightly  effected  hy  the  pressure  perturbation  but,  as 
already  known  previously,  are  greatly  influenced  by  the  entrainment  effect 
which  is  inversely  proportional  to  the  radius,  whereas  for  large  cloud 
the  reverse  is  true.  Titus  the  presence  of  the  optimal  cloud  size  when 
noth  of  these  two  factors  are  taken  into  consideration  appears  to  repre¬ 
sent  a  balance  between  these  two  opposing  effects. 

ilie  importance  of  the  initiaL  updraft  radius  and  velocity  is  shown 

in  Fig.  10  for  three  vertical  velocity  profiles  computed  using  a 

2 

constant  mass  flux,  M  =  u  R\  at  the  cloud  base.  The  three  profiles 

c  c 

differ  considerably  at  higher  levels  due  mainly  to  the  entrainment 
processes,  thus  indicating  especially  the  importance  of  the  proper  choice 
■'I  the  updraft  radius.  Squire  and  Turner  (1962)  in  contrast  did  not 
i 'Plain  this  large  variation  mainly  due  to  the i r  use  of  the  theoretical 
saturated  environmental  sounding.  For  radii  nr  ch  larger  than  those 
considered  here,  e.g.,  R  ■  10,000  m,  the.  importance  of  entrainment  is 
replaced  hy  the  importance  of  the  pressure  perturbation,  and  in  contrast 
i"  the  results  in  Fig.  10,  the  largest  radius  cloud  will  have  the  smallest 
vortical  velocity  maximum.  Nevertheless,  the  largest  cloud  will  have  the 
dm ■  pes t  peue l rat  ion. 

Manv  previous  investigations  have  attempted  to  account  for  the 
in l  I  lienee  of  the  pressure  perturbation  by  utilizing  a  reduced  gravity, 
i  e.,  ip  t lie  vertical  velocity  equation  g  is  replaced  by  g/(l+|-),  where 


.*  I 

if:  usual  Iv  assip.ned  .1  constant  value  ol  about  ().'>  and  t  lie  pressure 

perlurbat  i.011  term  is  otherwise  nop  lee  led  comp  I  o  t  e  I  v  As  pi  ey  i  011s  I  v 

pointed  out  tin-  pressure  pert  nrb.il  Ion  is  included  In  tills  model  since  1  lie 
vertical  variation  ol  0  is  retained  In  (lie  vert  Icilv  cquat  ion.  Hut  it 

In  Kqs  (.’Hi,  .’I*)  we  define 

•>  ) 

R‘  d  w 

b  -  ,  . 

a  w  d  /. 

tlu'vi  for  the  spec  i  a  I  case  if  equals  a  constant  t  lie  mod  I  I  i  ed  v.ravilv  loiuiu 

l.  it  ion  is  old  lined.  In  realltv  |>  Is  not  const  ant  as  shown  in  lie.  II. 

for  this  mean  tl  epical  sound  inp,  note  that  I  lie  reduct  ion  ol  the  oiled  ive 
Imovaticv  Is  ju'iie  ra  1  1  v  much  less  than  the  common  I  v  used  lactoi  .’/I, 

I'll’..  I.’  shows  anothei’  set  ol  |>  prollle;,  I  rom  an  Individual  Salem  III. 

sound  (np.  which  contains  several  stable  lovers  above  I  lie  cloud  base  at 
•i  h' .  <>•'(  mb.,,  as  indicated  bv  posit  Ive  values  el  AS  /,  1 /.  In  Kip..  I  t. 

I  lie  H  profiles  lor  K  .’000  and  ’>000  m  reprere.il  rd  bv  the  daslied  and  solid 
lines  in  Kip..  I  ,  t  espec  t  i  ve  I  v  ,  exhibit  tarp.e  variations  and  a  rep.  ion  ol 
nopal  Ive  values,  indicat  inp,  tli.it  the  pri'ssure  perturbation  lurcc  acts  to 
p.trtiallv  counteract  (he  buoyancy,  thus  smoothinp,  out  t  lie  veilic.il 
velocl.lv  pro!  ili";  as  shown  In  Kip.  1  a .  A  cloud  of  small  radius,  K  1000 

m.  Is  unable  to  break  t  hroup.h  the  stable  lovers  wbeieas  the  .‘000  and 
’'ODD  m  radius  clouds  are  capable  but  their  veil  leal  velocity  piotiles 
are  si  I  I  I  I  i  11 1  I  iienced  bv  it.  Willi  t  lie  pressure  perturbation,  as  re 
presented  bv  tin'  dashed  cmves,  the  i  input  t  ol  the  stable  I  avers  are  sub 

st  ant  i  .1 1  I V  reduced,  especially  lor  tin  >000  m  radius-  or  >0  ill  l.up.er  cloud 
The  perlurb.it  ion  pressure  f  1 1 1 1  ■  •  1  It  In  0111  solut  ions  can  be  solved  loi 


directly  by  integrating  the  traditional  vertical  velocity  equation  using 

our  computed  velocities.  The  vertical  profiles  of  the  scaled  perturbation 

pressure,  ir*  -  11  -  11  ^ ,  where  -n  is  the  ratio  of  the  perturbation  pressure 

divided  bv  environmental  densiLy,  is  found  to  be  negative  at  low  Levels 

and  positive  near  the  top  of  the  cloud  as  shown  in  Fig.  15  for  cloud 

radii  of  2000  (dotted  curve),  5000  (dashed  curve)  and  10,000  m  (solid 

line).  The  unknown  constant  of  integration  it  was  assigned  the  value  of 

o 

ii  at  400  mbs.  The  slope  of  each  curve,  i.e.,  Dir*/;tz,  indicates  the 
degree  to  which  the  buoyancy  is  adjusted  bv  the  non-hvdrostat  i  c  pressure 
perturbation.  Positive  slope  implies  that  the  buoyancy  is  reduced  whereas 
negative  slope  indicates  an  enhanced  buoyancy,  as  found,  for  example, 
near  Lhe  cloud  top.  The  Salem  111.,  n*  profile,  not  shown,  contains  a 
negative  slope  between  667.64  and  472.64  mbs,  showing  that  the  pressure 
perturbation,  in  part,  is  offsetting  the  reduction  in  buoyancy  caused  by  the 
stable  layers  found  in  the  sounding  between  7.32  and  572  mbs. 

In  Fig.  16  the  perturbation  pressure  profile  is  displayed  In  the 
same  format  as  m*  above.  Note  the  exceeding!'  small  vilues  of  less  than 
0.2  mbs  for  K  =  2000  m  (dotted  line)  and  the  large  positive  and  negative 
values  at  the  cloud  base  and  lop,  respectively,  for  a  radius  of  10,000 
m  (solid  line).  Also,  with  compress i ona 1  heating  (dashed-dot  curve), 
tlicri'  is  a  slight  reduction  when  compared  against  the  no  compress  Iona  1 
healing,  case,  as  represented  hv  the  dashed  curve. 

Vertical  profiles  ol  the  average  liquid  water  concentration  for 
c I ou.l  radii  ol  500  (Jolted  curve),  2000  (dashed  curve)  and  5000  m  (solid) 
line),  after  fallout,  are  shown  in  fig.  17.  It  is  seen  that  associated 
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with  increasing  cloud  size,  and  hence  larger  vertical  velocities,  is  (lie 
upward  movement  of  the  maximum  (j  toward  the  cloud  top.  A  similar  trend 
is  not  clearly  evident  in  the  fallout  rate  within  each  5  mb  interval  as 
illustrated  in  Fig.  18  for  the  same  cloud  radii.  Except  for  the  500  m 
cloud  (dotted  curve)  the  fallout  rate  is  a  maximum  at  nearly  the  middle 
point  of  the  vertical  profile.  The  total  precipitation  rate,  for  the 
tropical  sounding  without  ice  phase,  is  1.57.  4.73  and  5.17  cm  hr-1  for 
cloud  radii  of  500,  2000,  and  5000  m  respectively.  The  addition  of  tile 
ice  phase  will  increase  the  precipitation  rate  to  various  degrees,  rouchlv 
10-25  per  cent  depending  on  details  of  its  inclusion.  rti  these  cal¬ 


culations  the  constant  of  proportionally  A  in  the  fallout  term,  F,<) .  (7/), 

U  ”Ith  tot  100  .  <  R  <  2000  .  .  s«  10'V  for 

R  >  2000  m.  The  Salem  Ill.  sounding,  with  an  ice  phase  insetted  linearlv 
between  -10  and  -401  C  and  with  A  2  6.6  *  10  \  ,  produces  prec i p i t a t ion 

rates  of  0,69,  5.34  and  5.79  cm  hr  ^  for  cloud  radii  of  1000,  2000  and 
5000  m  respectively.  The  large  increase  between  1000  and  2000  m  is 
due  to  the  fact  that  for  the  larger  radius  the  cloud  is  able  to  peuotrale 
the  stable  layers  and  reach  a  considerable  height. 

It  should  be  pointed  out  that  the  precipitation  rates  given  above 
represent  the  average  precipitation  rites  for  the  developing  period  of 
the  different  types  of  clouds  only.  If,  instead  of  decaving,  the  deep 
cumulus  cloud  and  its  circulation  *-oinain  in  a  semi-steadv  state  after 
tully  developed,  then  the  preoi  pi  t  u  '  e'i  rate  can  be  much  higher  than  the 
values  given  above  since  then  almost  all  the  water  vapor  coming  into 


the  cloud  through  its  base  and  sides  will  be  precipitated. 
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6.  Influence  of  cumulus  convection  on  the  mean  properties  of  the 
at  mosphere . 

We  have;  illustrated  how  environmental  factors  near  the  updraft 

can  influence  the  basic  o Loud  properties.  The  reverse  process,  i.e., 

Lite  net  effect  of  short  lived  cumulus  convection  upon  the  thermodynamic, 

moisture  and  momentum  structure  of  meso  and  synoptic  scale  systems  has 

been  reeop.nl zed  for  years.  Various  parameterization  techniques  have 

been  introduced  to  incorporate  the.  vertical  distributions  of  condensation 

heating  and  momenLum  transfer  into  the  large  scale  flow,  see,  e.g., 

Kuo  (.196'),  1974),  Kreitzberg  and  Perkey  (1976)  and  Anthes  (1977). 

II  is  well  known  that  the  overall  effect  of  cumulus  convection  is 

to  warm  its  surroundings.  The  contribution  to  the  mean  heating  from  a 

single  cloud  is  from  the  combination  of  the  latent  heat  and  the  upward 

l ransporl  of  sensible  heat.  We  want  to  compute  the  mean  heating  rate, 

aT/'ii,  produced  by  llu>  cumulus  convection,  including  both  the  effects  of 

t he  compress  tonal  heating  in  the  far  environment  and  the  effect  of  the 

moist  and  cool  downdraft.  bet  us  represent  tin  areas  ol  the  updraft, 

moist  downdraft  and  far  away  downdraft  bv  A  ,  A,  ,.nl  A.,  respectively, 

c  dw  d 

and  let  the  vertical  velocities  In  these  areas  he  w  ,  w.  ,  and  w..  Then 

i  »  l W  ll 

tlu'  condition  of  zero  net  vertical  mass  flux  can  bo  taken  as 

Aw  b  A.  w  +  A  w  =  0  (10) 

c  c  dw  <lw  d  d 

whore  the  influence  of  the  difference  in  density  in  the  different  areas 

lias  been  disregarded  for  simplicity.  Further,  we  shall  express  w,  ,  " 

dw  uw 

in  tonus  ol  w  o'  and  A  of  the  cloud  updraft  bv  Introducing  the  non- 
c  c  c 

necal  ivo  factors  of  proportionality  u.,,  i  and  k  and  wi  iting  them  as 


and  \  , 
dw 
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w,  =  -a..w  ,0'  =  a0 '  ,  A  .  =  kA  . 

dw  2  c  dv  c  dw  c 


(  n  i.i..,  > 


On  substituting  these  relations  in  (50)  we  then  obtain 

Aw  =  -(1  -  u.k  )  A  w  .  (id ) 

d  d  2  c  c 

We  consider  that  the  heating  effect  of  the  cumulus  convection  in  (he 
layer  Az  is  produced  by  the  following  three  different  processes: 

i)  the  net  upward  transport  of  sensible  heat  into  this  layer  inside  1 
and  in  tin  convert ive  ,uo  ist  downdraft  in  the  immediate  environment 
i  i)  the  compress ional  heating  created  by  the  slow  descending  motion  ii 
the  far  away  environment; 

iii)  tire  release  of  latent  heat  in  the  cloud  updraft. 

Here  the  first  two  processes  are  directly  concerned  with  the  sensible  heat 


while  iii)  is  due  to  the  latent  heat.  We  shall  calculate  these  three  contri¬ 
butions  separately  below. 

According  to  the  relations  (31a, b,o)  the  area  average  of  the  upward 
heat  transport  T*  inside  the  cloud  and  its  immediate  moist  environment  is 


given  by 


where 


T*  -  ( I  -  u.jnk )  a  w  0', 
0  2  o  e  c 


ao  A 


(  11a) 


Thus,  on  using  R“  to  represent  A  for  the  ax isymmet r ic  cloud  we  find  that 
the  mean  heating,  effect  produced  by  the  first  process  is  gfven  bv 

"Vr’1  *=  * 

o 

The  contribution  ol  the  compress iona 1  heating  in  the  far  environment  t 
the  mean  heating  can  be  taken  as  given  by 

1,2  _d  d  _o  AT  , 

>t  An  Tz  '? 

n 
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when'  Af/t  is  the  ratio  of  the  existence  of  the  cloud  at  the  level  z 

A 


to  tl\e  total  time  t  of  existence  of  the  cloud.  Thus,  on  making  use  of  the 
relation  (12)  we  then  find 


:>  T . 


.  .  ..  a  w  „  a  [)  0 

1,2  _  , ,  o  c  At  o 

Vt  ~  (1  -  “2k)  —  i  rz" 

o 


(35) 


flu'  sum  of  these  two  effects  represents  the  total  heating  from  sensible  heat 
We  shall  represent  tills  sum  by  DT^/St,  viz.. 


.1  T 


:>  t 


pn 


aw  a  50 

/  ,  / 1  ,  ,  o  c  At  o 

(Pi*)  +  (1  -  «,,k) - r: 


Oz  0 


it  I  f)z 


(36) 


white  T*  is  given  by  (11) 


0 


The  heating  due  to  the  release  of  latent  heat  can  be  taken  as 
■  IT.,  .  3q 


I, 


*s 


,n 


"  •*o  c  wcl  i)—  "  "(q*  -  qs}1  ’ 
p 


(37) 


when-  q*  is  tiie  value  ol  q  just  outside  the  updraft.  from  the  results  given 


ihove  we  t  i nd  that  the  total  heating  rate  is  given  by 

)T,  ,)T.. 


IT 


a  t  a  t 


a  t 


(38) 


where  i  Ij/tt  is  defined  in  (36)  andilT^/Ht  in  (37). 


Nolo  Liu1  reduction  in  the  heating  rate  profile  shown  in  Fig.  19 
lot  a  1000  m  updraft,  computed  with  the  effect  of  a  moist  and  cool  down- 
dial  l  (dashed  1 i no )  as  opposed  to  that  computed  without  this  effect  (solid 
line).  Mere  the  ratio  of  moist  downdraft  area  to  the  updraft  area. 


i.e.,  k  A,  /A  ,  is  taken  as  1.0,  and  the  ratio  of  vertical  velocities, 
dw  c 


-w  /w  ,  is  taken  as  O.A.  In  both  cast's  the  role  of  compress  iona  1 
dw  c 


heal  ini'  has  been  included,  assumming  that  the  ratio  ot  updralt  area  to 

lot  il  ana,  i.e.,  A  /A  a  is  0.1.  The  total  sensible  heating  profile 
co 

lor  l  lie  moist  downdraft  case  ( dashed -do t  curve),  is  slightly  negative  at 
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low  levels  and  only  slightly  positive  at  high  levels.  Thus  latent  pro¬ 
cesses  are  responsible  for  most  of  the  total  heating.  The  contribution 
of  the  corapressiona 1  heating  in  the  environment  to  Lite  total  sensible  heat 
is  indicated  by  the  dotted  curve. 

The  significance  of  latent  heating  (dotted  curve)  is  shown  in  fig. 

20  for  a  2000  m  radius  cloud  containing  compressional  heating,  along 
with  the  total  sensible  heating  rate  (dashed  curve)  and  their  sum,  i.e., 
tiie  total  heating  rate  (solid  line).  Note  the  importance  of  the  total 
sensible  heat  contribution  near  the  cloud  top  but  otherwise  the  heating 
is  due  mainly  to  latent  heat.  This  peaking  near  the  cloud  top  is 
even  more  pronounced  In  the  'cloud'  sensible  heating  rate  (dashed-dot 
curve.  Fig.  21)  for  a  5000  m  radius  cloud.  The  magnitude  of  the  negative 
contribution  at  lower  levels  is  also  fairly  large  but  is  considerably 
diminished  provided  the  positive  contribution  from  compressional  beating 
(dotted  curve)  is  added  to  form  the  total  sensible  heating  rate.  The 
importance  of  compressional  heating  is  revealed  also  in  Fig.  21  by  the 
difference  between  the  total  heating  rate  computed  with  environmental 
effects,  a  =  0.4,  b  =  -c.  /I.it,  (dashed  curve)  as  compared  wiLh  the  'cloud' 
total  heating  rate  (sol  ill  line).  Note  that  this  difference  is  smaller 
in  magnitude  than  the  contribution  from  compressional.  heating  itself 
(dotted  line).  The  reason  why  is  because  compressional  heating,  also  tv- 
dui es  the  vertical  velocity,  thereby  reducing  the  local  rate  of  latent 
heat  released.  This  effect  is  Illustrated  further  in  Fig.  22  for  three 
different  cloud  radii.  Some  measure  of  t  lie  increase  of  the  total  heat  ing  rate 
with  increasing  cloud  radius  can  also  he  gained  by  comparisons  between  the  solid 


For  clouds 


id 


nirvi'S  in  Figs.  19  and  20  and  the  dashed  curve  in  Fig.  21. 
of  radius  R  -  1000  m  the  total  heating  rate  maximum  is  located  nearly 
midway  in  the  vertical  profile,  huL  it  is  shifted  upward  toward  the 
cloud  Lop  as  the  radius  is  increased. 

We  have  shown  that  tor  large  cloud  radii,  e.g.,  K  '>000  m,  com¬ 
press  tonal  heating  increases  the  total  heating  rate.  This  conclusion  is  also 
valid  Cor  cloud  radii  of  2000  and  500  m  (not  shown),  hut  obviously  pro¬ 
portionally  less.  For  the  range  of  parameters  considered,  the  role  ol 
the  moist  and  cool  downdraft  is  found  to  increase  the  sensible  heat  trans¬ 
port  T*  (Eq.  33)  slightly  by  the  addition  of  the  small  positive  lactor 
(-a.^ak)  in  T*  but  also  shown  to  reduee  the  vertieal  velocity,  so  that  there 
is  some  net  reduction  in  the  total  heating  rate  by  this  el  loot. 

It  should  be  pointed  out  that  the  expression  (37)  lor  the  heating, 
rate  from  latent  heat  is  based  on  the  total  condensation  rati'.  When 
the  cloud  dissipates,  the  Liquid  water  carried  by  the  cloud  will  re¬ 
evaporate  and  therefore,  the  mean  net  heating  rate  produced  by  the  cloud 
during  its  whole  life  span  should  be  given  by 


(e) 

1T  3T  I  T  A  l,q  ’ 
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where  q^^  stands  for  the  amount  of  liquid  water  re-evaporatei 
l  w 
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The  rate  of  moistening  of  the  atmosphere  is  given 


Ac  ,  .  (e), 

JC  =  A  (qs  -  ‘<o  +  qiw  } 


Because  of  our  lack  of  more  reliable  knowledge  of  the  relat  ion 


(  19) 

into  the 
hv 

(10) 

between 


the  cloud  tvpe  and  ,  the  intluence  in  re-evaporation  ol  the  cloud 
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water  into  the  environment  has  not  been  included  in  our  ca I  in  1  at  ions  vot. 

It  we  assume  that  all  the  unprec ip i t a t rd  liquid  water  drops  in  the  cloud  .  i. 
re-evaporated  into  tin  environment  at  the  same  level  as  thev  ooeur  in  the 
eioud  as  the  eloud  dissipates,  then  a  larger  cooling  will  take  plane  near 
the  eloud  top  than  at  low  levels  as  'lie  eland  dissipates. 

It  is  evident  that  the  various  area  mean  clients  produced  hv  a  sins  I! i. 
tvpe  of  cumulus  can  he  calculated  by  the  use  of  the  cloud  ph  .'s  I  pi  ,yi  led 
the  area  ratio  a  ^  (=A  /A  )  .ind  the  other  phvsica!  parameters  are  known. 
Obvious  Lv,  a  model  cloud  representative  of  Ihe  average  c, audit  sms  can  aka 
be  constructed,  by  combining  a  number  of  clouds  of  different  radii.  I'lie 
heating  and  transporting  influences  of  such  a  cloud  can  !>c  obtained  linn 
those  ot  the  individual  clouds  by  assigning  proper  weight  ini’  factors  to  then; 

The  use  of  the  mean  stability  and  boundary  laver  propel  ties  oi  the 
itmosphere  and  large-  and  meso- scale  lifting  to  determine  whether  cumulus 
convection  will  occur  or  not  in  the  region,  and  to  use  the  large-  and  mesa- 
scale  net  convergence  of  moisture  as  the  measure  of  the  average  candeiisat  Ian 
rate  in  such  systems  as  proposed  hv  the  first  author  (Kuo,  I  "■><>  > ,  1  *l  7  •'*  > 
still  appear  to  be  the  proper  road  for  the  paramatic  t  epresent  at  ion  el  the 
time  and  space  average  heating  effect  of  cumulus  convection.  However,  the 
question  of  whether  cumulus  Convection  will  occur  at  a  given  locality  and 
a  given  time  can  not  be  answered  hv  utilizing  large  scale  propc  rt  ion  alone. 
The  problem  is  complicated  further  bv  the  fact  that,  cumulus  convert  ion  is 
of  the  nature  of  overturning  and  lienee,  it  restricted  to  a  small  region,  it 
can  occur  only  once  and  after  completing  this  process  the  stratification 
will  become  stable  and  the  moisture  content  exhausted.  On  the  other  hand. 


It 

f 

F 

t 


i!  the  convection  lakes  an  organized  form  with  a  cont  inuous  supply  ol 
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moist  .lit-  i  ruin  I  ar  ,iwav,  sucii  as  t  host*  in  the  muso-sea  1 e  storms,  then 
tin  system  eau  have  a  nineli  longer  1  i  t  e ,  Kvulentlv  whet  lie  r  such  a  system 
Kill  he  created  in  a  given  region  or  not  can  only  he  determined  bv  studying 
tin  aaierai  .Ivnamie  and  t  he  rnsodynam  i  c  conditions,  and  not  hv  using  t  lie  mean 
l.uve  stall1  st  rat  i  l‘ ica  t  ion  and  transport  properties  alone. 

In  view  ot  the  mans'  uncertainties  concerning,  the  average  heating  effort 
>!  *  mini  Ins  convert  Lon ,  it  appears  preferable  to  determine  it  bv  a  less  i n— 
ml  veil  simple  method  rather  than  by  the  use  of  a  detailed  model.  One  simple 
ve,  is  to  use  tiie  vertical  distribution  ol  the  potential  temperature  excess 
’  •  i  ceil  bv  I  1 S )  or  !'  given  by  (I8h)  as  the  vertical  distribu!  ion  ot  the 

heal  ins  rite.  on  ■  omp.  i  r  i  ii  g  the  total  healing,  profiles  in  l’icures  14,  .it), 

•  .,)  ’1  i or  clouds  with  radius  greate r  than  i 000  m  with  the  eloud  potent  ia! 
ic  I'liatme  pro  tile  in  fig.  l>  we  see  that  l  hey  are  lairlv  similar  to  each 
M,i'  r.  indie  it  Lug  t  hat  the  d  i  st  r  i  hut  ion  ot  the  potential  temperature  excess 
'  in  be  used  as  a  reasonable  approximation  to  the  vet  t  ica I  distribution 
,  :  :  in  mean  heating  el  tort  created  by  the  eloud  aetivitv,  .e  proposed  bv 

■  '  ! .  !  I  ’ '  I  r  l  )  . 

In  t }  Hi.  a  ion  ol  how  tlu  latent  neat  ..id  <  ns  i  h  1  >  b.e.it  accunu 
0  u  all  individual  .  loud  is  transiiutted  to  tin  mean  atmosphere  nav  1  so 
•••  Ml. :  I  ■  C.c.  I .  Hot  i-  ot  least  two  dill  oreilt  pi  oeesses  can  i>e  visual  i?.ed. 
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7.  Summary. 

A  quasi-one-dimens ional  and  quns i-Lagrang ian  steady  state  modi’ I 
which  takes  into  consideration  the  influences  of  the  pressure  perturba¬ 
tion,  vertical  density  stratification  and  nearby  moist  downdraft  on  i  lie 
axi symmetric  cumulus  cloud  convection  is  obtained  by  integrating,  the 
combined  heat  and  moisture  equations  and  the  azimuthal  vorticltv  equation 
from  the  cloud  center  to  the  outer  edge  and  utilizing  verfous  entrain¬ 
ment  considerations.  This  model  shows  that  the  Influence  of  the  pressure 
perturbation  is  proportional  to  the  square  of  the  cloud  radius  and  henci 
it  is  small  for  small  clouds.  This  model  is  used  to  investigate  the 
variations  of  tiie  various  cloud  properties  with  the  radius  of  cloud 
under  two  different  environmental  conditions,  one  is  for  the  potent  i  a  1 1 v 
more  unstable  and  more  smooth  mean  tropical  atmosphere  during,  the 
hurricane  season  while  the  second  is  for  a  less  unstable  middle  latitude 
environment  with  prominent  thin  stable  1  avers  in  mid-lower  troposphere, 
further,  two  different  kinds  of  cumulus  convection  are  considered, 
namelv.  A)  is  of  the  form  of  a  simple  updraft  surrounded  by  the  almost 
unaltered  iar  environment,  either  with  or  without  a  feed  back  from 
the  compt  essional  boating,  in  this  environment,  while  I?)  is  with  a  nearly 
saturated  and  cooler  nearby  downdraft,  as  is  often  observed  in  the  cloud 
clusters  In  the  tropics.  In  tliis  latter  case  the  cloud  is  not  influenced 
by  the  far  environment  direct iv. 

It  Is  found  that  t ho  pressure  perturbation  in  the  cloud  can  roach 
from  +  2,5  to  +  7  mb  at  the  cloud  base  ami  the  cloud  top  as  the  radius 
of  tilt  cloud  i ti»- reuses  from  5  km  to  1(1  km  but  is  small  for  cloud  radius 
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less  than  2  km.  Further,  the  magnitude  of  the  perturbation  pressure  is 
reduced  hv  the  Inclusion  of  the  temperature  feed  back  from  the  compres- 
sionallv  lieated  environment. 

Since  the  temperature  equation  (18)  involves  neither  the  pressure 

nor  the  vertical  velocity,  the  pressure  perturbation  has  no  direct 

influence  on  Lhe  cloud  temperature  excess.  On  the  other  hand,  the 

inclusion  of  the  pressure  perturbation  reduces  the  vertical  velocity  in 

Lhe  cloud,  as  it  must  since  pressure  perturbation  transmits  energy  to 

the  environment.  One  interesting  and  possibly  also  important  effect  of 

the  pressure  perturbation  is  that  it  gives  rise  to  an  optimal  cloud  radius 

of  about  8  km,  corresponding  to  the  highest  value  of  w  ,  under  the 

max 

tropical  environment,  while  without  the  pressure  effect  w  increases 

max 

mount  on  Lea  1 1 v  with  increasing  R. 

The  influence  of  Lite  density  stratification  is  to  reduce  w  slightly 

below  the  level  of  w  and  to  increase  w  above  this  level  and  also  to 

max  c 

increase  the  cloud  depth  slightlv,  but  the  total  influence  is  not 
prominent.  The  results  show  that,  for  the  norm  i  tropic  il  atmosphere, 
lhe  cloud  properties  change  gradually  with  the  cloud  radius,  while  for 
the  less  unstable  environment  with  prominent  stable  layers  in  mid-lower 
stratosphere,  I  lie  cloud  properties  change  abruptly  as  the  cloud  radius 
increases  t tom  loss  than  to  larger  than  J  km,  reflecting  the  inhibiting 
ini  I  nonce  of  the  shallow  stable  I  aver  in  mid-lower  troposphere'  on  the 
‘oilier  clouds. 

It  i  s  also  tomui  that  the  presence  of  a  moist  hut  eoo  ler  cotivec  t  i  ve 
wn  !  t  ilt  i educt-s  Pot  li  t  he  updraft  and  the  cloud  depth.  It  also 
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augments  the  prec i pi t at  ion  rate,  especially  Cor  the  smaller  clouds. 

The  heating  and  moistening  influences  of  cumulus  cloud  convection 
on  the  large  scale  atmosphere  have  also  been  calculated  by  the  model.  It 
was  found  that  the  main  part  of  the  heating  is  from  the  latent  heat  but 
compressional  heating  of  the  environmental  air  is  also  of  importance, 
especially  close  to  the  top  of  the  cloud.  On  comparing  the  average  heat  in 
rate  with  the  potential  temperature  excess  6 '  given  bv  this  model 
we  find  a  close  similarity  between  them  for  cloud  radius  larger  than  1  km. 
It  is  therefore  concluded  that  3'  can  be  used  as  a  rough  approximation  to 
the  vertical  distribution  of  the  average  net  heating  effect  produced  hv 
cumulus  convection  in  the  parameterization  scheme.  Of  course,  results 
obtained  from  detailed  calculations  from  the  cloud  model  can  also  he  used. 

An  equation  for  the  vertical  mixing  of  the  horizontal  momentum  by 
cumulus  convection  on  the  change  of  the  horizontal  velocity  has  also  been 
obtained . 
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Appendix:  List  of  Symbols 

)  radial  distance  from  the  axis  and  altitude 

radial  and  vertical  components  of  velocity 

mean  vertical  velocity  of  cloud  updraft 

,  ,  T  ,  0  ,  (\  :  pressure,  density,  absolute  temperature,  potential 

o  o  o  o 

temperature  and  water  vapor  mixing  ratio  in  the  undisturbed 
environment,  all  taken  as  functions  of  the  altitude  z  only 
T' ,  0*  (=  nT' ) ,  <| ' :  mean  values  of  departures  of  the  total  pressure 
p,  temperature  T,  potential  temperature  0  and  water  vapor 

mixing  ratio  q  in  the  cloud  from  the  undisturbed  values  p  , 

o 

T  ,  0  and  <i  . 
o  !  >  o 

-  T ( 1  +  O.bOhq):  virtual  temperature  of  moist  air 
■  ii T  virtual  potential  temperature 

specific  heat  at  constant  pressure 

1)  gas  constant  of  dry  air;  2)  radius  of  cloud  updraft 

1000  mb,  R/c  _ 

\  -  )  P  i 


l^ti),  (|  saturation  mixing  ratio  at  temperatare 

1 1'inpe  rat  u re  T  of  environment 
o 


<l”  -  <10(T)-(| 


cloud  and  at 


q  <1,  +q  I  iiuiid  w.iLci  mixing  ratio  in  cloud  updraft;  q  is  cloud 

.  '  C  W  t  W  1  v‘\\' 

liquid  water  mixing  ratio  and  o  is  mixinc,  ratio  of  rain  water 
1  i  w 

i ’  rat  e  o I  condensa t  i i >n 

I  latent  lie.it  of  condensation,  which  mav  also  include  the 


latent  heat  of  freezing 
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Fig.  1.  Profiles  o  I  environmental  potential  temperature  0  (solid  line) 
and  equivalent  potential  temperature  0  (dashed  and  dotted  lines) 
for  the  mean  Culf  of  Mexico  hurricane  season  sounding.  In  the 
ilotted  curve,  variation  of  latent  heat  of  vaporization  with  tem¬ 
perature  is  allowed  and  latent  heat  of  freezing  is  included  linearly 
between  -10°  anil  -A0*C, 

Fig,.  ’ ,  Vertical  profile  o  f  )  S  /i!z  for  the  tropical  sounding  in  Fig.  1. 

Fig,,  i.  Cloud  temperatures  excesses  for  radius  R  1000  m  witli  (a  =  -0.  I  ,  b-0. 00V , 

dashed  curve)  and  without  (a  =  0.0,  solid  line)  moist  downdraft  cooler 
than  the  undistributed  air. 

Fig,.  A .  Vertical  velocity  profiles  for  cloud  radius  R  "  1000  m  with 
tdashed  curve!  and  without  (solid  line)  moist  downdraft.  The 
dotted  curve  includes  tlu‘  small  add  i  t  ioilal  influence  of  evaporation 
at  the  cloud  houndary  (a  ,  -•  0.1  r>)  for  the  downdraft  case. 

Fig,.  Same  e  I  end  as  Fig,.  A  except  showing,  t  I"-  add  i  t  ienal  influence  el 

the  latent  heat  of  fusion  upon  t ho  vertical  velocity.  Dashed  curve 
includes  latent  heat  of  freezing  whereas  the  solid  curve  does  not  . 

F  ig .  <i.  Cloud  potential  temperature  excesses  for  three  cloud  radii  with  t.  a  0 .  >  , 

l>  (•  '  I .  :i  ,  dashed  curve!  and  without  comp  ross  i  onn  I  heal  iiir  (a  t'.O, 

P 

solid  line!.  No  nearby  moist  downdraft. 

I  ig.  i .  Vert  foal  voloeit  i  es  for  a  r>00  m  radius  cloud  computed  with  (a  t  > .  •  * .  h 
|> 

,  dolled  curve!  and  without  (a  0.0,  solid  line!  compress i ona I 
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heal  ills  and  pressure  perturb  it  ion  el  tool  s  (dashed  curve!  . 
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Fig.  8.  Same  as  Fig.  7  except  for  a  5000  m  radius  cloud.  Mash-dot t ed 
curve  is  without  compress ional  heating  but  with  both  the  o  term 
in  (21a)  and  pressure  effect  included  and  is  to  hi1  compared  witli 
dashed  curve. 

Fig.  9.  Variation  of  w  with  cloud  radius  with  (dashed  curve)  and  with- 
p  max 

out  (solid  line)  pressure  perturbation. 

Fig,  10,  Vertical  velocity  profiles  for  3  different  cloud  radii  taking 

the  cloud  base  mass  flux  as  constant.  1)  w  =  2  ms  '  and  R  =  2000  (2) 

o 

(dotted  curve);  2)  w  =  1  ms  ^  and  R  =  2000  m  (solid  line);  1) 

o 

-1  k 

w  =0.5  ms  and  R  =  2000  x  2 ~m  (dashed  curve), 
o 

Fig.  11.  Vertical  profiles  of  gravity  reduction  factor  d  for  R  300  m 
(dotted  curve),  R  =  2000  m  (dashed  curve)  and  R  =  5000  m  (solid 
line) . 

Fig.  12.  Same  as  Fig.  11  except  for  Salem  Ill.  sounding  in  Fig.  i  1. 

Fig.  13.  Vertical  profile  of  DS  /dz  for  a  1  BOO  C.S.T.  Salem  111.  sounding 

os 

(8/22/75)  . 

Fig.  14.  Vertical  velocity  profiles  for  3  radii  obtained  with  the  Salem, 
III.  sounding.  Dashed  curves  include  pressure  perturbation,  solid 
lines:  without  pressure  effect. 

Fig.  15.  Variation  of  perturbat ion  pressure  to  density  ratio  n*  with 

height  in  cloud  for  mean  hurricane  season  sounding,  for  cloud  radii 
R  -  2000  (dotted  curve),  K  -  3000  (dashed  curve)  and  R  10,000  m 
(solid  line).  No  compress  Iona  1  heating. 

Fig.  lb.  Profiles  of  the  perturbation  pressure  plotted  in  tin-  same 

format  .as  a*  in  Fig.  15,  will)  the  add  it  ion  ot  the  dashed-dot  eurve 
computed  including  compress  ional  heating,  for  a  5000  m  cloud. 
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Kin.  17.  Vortical  profiles  of  the  liquid  water  concentration,  q  ,  for 

cloud  radii  of  500  (dotted),  2000  (dashed)  and  5000  m  (solid  curve). 

Fig.  to.  Same  as  Fig.  17  except  showing  the  fallout  rate  within  each  5 
mb  interval. 

l'ig.  19.  Heating  rates  for  cloud  of  radius  1000  m,  a  =  0.1  k  =  1.0, 

o 

and  u,  =  0.4,  computed  with  (a  =  -  0 . 1  , b- 0. 002 ,  dashed  curve)  and  without 
(a  =  0,  solid)  moist  downdraft.  The  total  sensible  heating  profile 
is  indicated  by  the  dashed-dot  curve  and  the  contribution  (if  com¬ 
press  ional  heating  by  the  dotted  curve. 

Fig.  20.  Vertical  profiles  of  the  total  sensible  heating  rate  (dashed 
curve),  the  latent  heating  rate  (dotted  curve)  and  their  sum 

o 

(solid  line)  for  a  2000  m  cloud  with  compress ional  heating  (a  =  0.4  h=-  ~) 

Fig.  21.  Heating  rate  of  a  5000  m  cloud.  Dashed  curve:  total  heating 

rate  with  compress  i  nna  1  heating  in  environment  and  n~0.4,  b=-e  /Lit  : 
solid  curve:  total  heating  rate  produced  by  cloud  alone;  Dash- 
dotted  curve:  'cloud'  sensible  heating  under  influence  of  com- 
pressional  heating  in  environment;  Dotted  curve:  contribution 
from  compress ional  heating  in  environment  alone. 

Fig..  22.  u  times  the  local  rate  of  latent  heating  for  1  different  cloud  radii  with 
(a  -  0.4,  I>  /lut  dashed  curve)  and  without  (a  =  0.0  solid  line) 

compress i ima 1  heating,  withouL  moist  downdraft. 
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